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Abstract 

The creation of massless scalar particles from the quantum vacuum by spherical 
shell with time varying radius is studied. In the general case of motion the equations 
are derived for the instantaneous basis expansion coefficients. The examples are 
considered when the mean number of particles can be explicitly evaluated in the 
adiabatic approximation. 
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1 Introduction 



The Casimir effect is one of the most interesting manifestations of nontrivial properties of 
the vacuum state in a quantum field theory (for reviews see [1], 0, ||, |J) and can be viewed 
as a polarization of vacuum by boundary conditions. A new phenomenon, a quantum 
creation of particles (the dynamical Casimir effect) occurs when the geometry of the sys- 
tem varies in time. In two dimensional spacetime and for conformally invariant fields the 
problem with dynamical boundaries can be mapped to the corresponding static problem 
and hence allows a complete study (see 0, |J and references therein) . In higher dimensions 
the problem is much more complicated and is solved for some simple geometries. The vac- 
uum stress induced by uniform acceleration of a perfectly reflecting plane is considered in 
|5|] . The corresponding problem for a sphere expanding in the four-dimensional spacetime 
with constant acceleration is investigated by Frolov and Serebriany || 0] in the perfectly 
reflecting case and by Frolov and Singh |§ for semi-transparent boundaries. For more 
general cases of motion by vibrating cavities the problem of particle and energy creation 
is considered on the base of various perturbation methods || |TTJ, [12], [13|, [L4], [15], |16| (for 
more complete list of references see ]16[). It have been shown that a gradual accumulation 
of small changes in the quantum state of the field could result in a significant observable 
effect. A new application of the dynamical Casimir effect has recently appeared in con- 
nection with the suggestion by Schwinger [17] that the photon production associated with 
changes in the quantum electrodynamic vacuum state arising from a collapsing dielectric 
bubble could be relevant for sonoluminescence (the phenomenon of light emission by a 
sound-driven gas bubble in a fluid [@). For the further developments and discussions this 
quantum- vacuum approach see [^1], |19|, [22], [23|] and references therein. In the present 
paper we consider particle creation from the quantum scalar vacuum by expanding or 
contracting spherical shell with Dirichlet boundary conditions. In next section we derive 
the equations for the instantaneous basis expansion coefficients and the formula for the 
Bogoliubov coefficients. The examples when this coefficients can be explicitly found at 
adiabatic approximation are discussed in section 3. 



2 General consideration 

Consider a scalar field p satisfying Dirichlet boundary condition on the surface of a sphere 
with time-dependent radius a = a(t): 

d 2 

(—-A)<p(x,t)=0, <p\ r=a{t ) = 0. (1) 

The corresponding eigenfunctions can be expanded in a series with respect to the instan- 
taneous basis 

<Pa( x ) = Qap(t)(j>p(x, a(t)), (2) 

13 

where the collective index a denotes the set of quantum numbers specifying the solution. 
Here 0/j(x, a) exp(— iuot) are the corresponding eigenfunctions for a static sphere with 
radius a: 

<Pl3(x,a) = 3/2-//- \ Y im(0,<P), (3) 

a Jl{Jl,n) 

(3 = (l,m,n), i = 0,1,2,..., —l<m<l, n = l,2,..., 
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with j\ n being the n-th zero for the spherical Bessel function ji(z), ji(ji, n ) — 0, Yi m (9, <p) 
is the spherical harmonic. Note that the sets of quantum numbers are the same for static 
and dynamic cases. Function (Q) satisfies boundary condition ([I]) . Putting expression (|2|) 
into the field equation (p]) we obtain 

Yl {foiQapit) + Up(t)q a p(t)] + 2^0/3 + q a p4>p} = 0, u fi (t) = ji, n /a(t), (4) 

P 

where dot stands for the time-derivative. Let us multiply this equation by (j)% and in- 
tegrate over the region inside a sphere at a given moment t. Using the orthonormality 
relation 

J <f>p<l>*pcPx = 8pp, (5) 

this yields 

q a p(t) + u 2 p (t)q a p(t) = J2(2q Ql3 >gp>(3 + q a p>g%), (6) 

P' 

where we have introduced notations 

<,*j3„, „(1) / A. A* 



gp'P = - J <t>p>4>*pd x, gyp = - J <pp<p* d x, (7) 

with integrations over the region inside the sphere. The relation between coefficients (|7|) 
can be found by making use the completeness condition for eigenfunction (|3|): 

J2M^(r) = $(r-r)- (8) 

7 

This yields 

(l) dgpp „ (n\ 

m = —^ + 2^9P'i9pr ( 9 ) 

By taking into account this relation we note that the general structure of equations (||) 
for the instantaneous basis expansion coefficients is similar to that for the plane case. 
Using the standard orthonormality relations for the spherical harmonic it can be easily 
seen that 

gp'P = g l n>n$ll'Smm', gpp = Srf^l^mm' • (10) 

To evaluate the integrals in (|7j) we use the formula for the integrals involving the Bessel 
spherical function 



aj'i(a)ji(b) - bii{a)][{b) 



z 2 3l (az) 3l (bz)dz = (ii) 

cr — cr 

From this formula the following expression for the coefficients can be obtained 



l _ i l (1)1 _ i l | 12 V 1 l l h — — fi9\ 

p a 

where 

<L =0, (13) 



'2 '2 ' 

3Lti' 3l,n 
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As we see the coefficients gp>p> are antisymmetric. This is a direct consequence of orthonor- 
mality relation (|5[). From expressions (^) it follows that the coefficients q a p in (^|) can be 
chosen diagonal with respect to the quantum numbers I and m, and are independent on 
m: 

Qpp' — Su'S mm 'q l nn ,. (14) 
Now the expansion for the eigenf unctions takes the form 

oo 

<Plmn = £ im><l>lmn>{x , «(*)) , (15) 
n'=l 

with the following infinite set of coupled differential equations for the coefficients: 

oo oo oo 

Qnn' + ^In'i^Qnn' = 2 ^ £ Qn P a pn' + ^ £ 1n V a V n> + ^ £ Qn P a P s a n's- (16) 

p=l p=l P>s=i 

If the sphere is asymptotically static at past and future then the in- and out- vacuum 
states can be defined by using the solutions for coefficients corresponding to the in- and 
out-modes cp^(t), ^fa Ut K^) with asymptotics 

- e'^5 a p, t^-oc (17) 
&\t) - e^'Sap, f-oo, (18) 

where we use the notations 

o% = —, < Ut = — , o± = lim ait) (19) 

for the corresponding eigenfrequencies. The field operator in the Heisenberg representa- 
tion may be expanded in terms of the corresponding eigenfunctions 

= E9$(*)M*> <*(*))> * = in ^t (20) 

P 

as 

v(^*)=E[°?V? ) + a? )+ ri' ) 1- (21) 
The in and out vacuum states |in >, |out > are defined in accordance with 

a$\s>=0. (22) 
The corresponding eigenfunctions are related by the Bogoliubov transformation 

# = EM° Ut) +^° Ut) *)> (23) 
P 

with Bogoliubov coefficients a a p and j3 a p- Substituting instantaneous basis expansion, 
multiplying by ip* and integrating over the region inside the sphere we obtain the corre- 
sponding relation between expansion coefficients: 

4 n) w = £K4; ut) + MF'*)- (2 4 ) 
p 
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Here we used the formula 

L* 1* j3 



4>p(j)*p,d X = 5nn>bll>bm-m>, (25) 

and the fact that expansion coefficients are independent on m. Now by taking into account 
relation (|lj|) the Bogoliubov transformation can be written in the form 



&\t) = u<$ t)l + pU P T )l y (26) 
p 



In particular, taking into account ([T7]), in the limit t — > oo from (^) we receive 

<(*) = W*^* + (3 aP e^r\ t +oo. (27) 

The mean number of out particles with quantum number a in the in-vacuum state is 
determined by the Bogoliubov coefficient (3 aa as 

< in|i\yin >= £ l/^l 2 , /3 m = (28) 

cr 

with the standard Klein-Gordon scalar product (see, for instance, |4]]). From (|28|) it follows 
that /3 a ' Q = P l n , n 8ii>5 m - m i and hence the total number of created scalar particles is given 
by 

oo oo 

<in|iV|in>=£ <in|iV imri |in>=^(2/ + l) ]T \(3 l n , n \ 2 (29) 

Imn 1=0 n,n'=l 

The coefficients a l np and in (^|) are determined from the solutions of infinite set of 
coupled differential equations ([16]) with time dependent coefficients. The problem can be 
simplified in limiting cases when various approximations can be used (see, for example, 
@-[|nj for the case of plane boundaries). In next section we consider the adiabatic 



approximation. 



3 Adiabatic approximation 

From the form of equations (|T^) it follows that there are two types of effects which lead to 
the particle creation (see also |15|]). The first one, called squeezing of the vacuum, is due 
to the nonstationary eigenfrequencies u)i n (t) as a result of a dynamical change of the radius 
of the sphere and is described by the second term on the left of (|I6|). The second one, 
referred as acceleration effect, is due to the motion of the boundary and comes from the 
terms on the right of flTB]) . Due to the antisymmetry of the coefficients a l nn , the squeezing- 
and acceleration- effects give additive contributions to the number of created particles per 
mode to first non-vanishing order of perturbation theory. In this section we will consider 
the squeezing contribution to the number of particles. Note that when the sphere is 
moving on a time scale slow on the scale of the created quanta frequencies (adiabatic 
approximation) (see |l ij and references therein) the terms on the right of ([IB] ) containing 



the derivatives of the sphere radius are small and the squeezing effect is dominant. In the 
adiabatic approximation the matrix q a p can be chosen in diagonal form q a $ = q a S a p and 
Inn 1 = Qin$nn' and from fll6|) we receive 

q ln (t)+uf n (t)q ln (t) = 0. (30) 
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As it follows from ( pjl) the necessary conditions for the adiabatic approximation are 

d 2 , ad<j 2 n . (31) 
As an example let us consider an exactly solvable case when 

a(t) = —== 1 A>\B\, (32) 

IA + Btanh(t/t ) 



where A, B and to are constants. This motion corresponds to the sphere contraction for 
B > and expansion for B < 0. The corresponding frequencies are 



Win(t) = x„VA + Stanh(t/to), (33) 

< = -j^=Ji, n VA^B, = -j^— = j ln y/A^B. (34) 

a — oo o oo 



Now we need to solve the equation ( j30|) with u>i in (t) given by fl33|) . The corresponding 
solutions are given by hypergeometric function. The normalized in- and out- modes are 
given by formula H 



qf n {t) = (2^ n )- 1 / 2 exp[-^+t-^oln[2cosh(t/t )]] 



x 



x 2 F 1 (1 + z^ n t ,^ n to;l=F^/n t o;^(l±tanh(t/t ))), s = in, out, (35) 
where uper/lower sign corresponds to the in/out- modes, and 

^ = ^KT±^n)- (36) 

Using the diagonality for the coefficients q l nn , over nn' we conclude that the same is the 
case for Bogoliubov coefficients in fl26|) : 



a nri — a ln^nn', Pn n ' — An^nn'- (37) 

Now the Bogoliubov transformation (00) can be written as 



„M _ _(out) . a _(out)* /oo\ 

Using the linear relation between hypergeometic functions, similar to [|J for the coefficients 
in this formula one finds 

„ ^rV /2 r(i -^to)r(-^rto) 
/n UlnJ r(-^o)r(i-^o))' 1 j 

, _ /< ut y /2 r(i-^t )rK°n%) {m 

Pln [ufcj ^(^u^ n tom + ^u J ^ n to)^ 1 ] 

The mean number of particles produced through the modulation of the single scalar mode 



is 



, sinh 2 (7ru; Zn toj 



< in|jV Iri |in >= \f3 ln \ 2 = ■ w in . ; ■ It out+ y (41) 

smh(7ru;^to) smh(7u^%) 



The total number of particles produced is obtained by taking the sum over all the oscil- 
lation modes inside the sphere: 



< in|tf|m >= ±(21 + 1) I ^iM^-VA^Bm 



(42) 



Therefore the energy related to the particles production inside the sphere is given by 

smh 2 [nj l>n t (V^+B ~ VA^B)/2] 



1=0 n=l 



i=o n=i swh(irj l>n t y/A - B) sinh(7r/' /in t v'^ + B] 



Ji 



, n VA + B. 

(43) 

It can be seen that the similar results may be obtained for another type of boundary 
conditions, for instance, for the Neumann or more general Robin conditions taking instead 
of ji tTl the corresponding eigenvalues. 

As an another example let us consider the case in which the sphere radius varies as 
(for the corresponding problem in the case of parallel plates see WW ) 



a(t) = a 



1 + 



cosh (t/tc 



-1/2 



(44) 



with a , b, t being constants. This function describes a pulse with characteristic duration 
to and modulation strength b. As in the previous example the corresponding in- and out- 
solutions to equation fl3"D| ) can be expressed in terms of hypergeometric function (see, for 
instance, [0]). Expanding in-modes in terms of the out functions for the mean number 
of scalar particles with a given mode one receives 



< in|iV; n |in > 



cos 2 [(7r/2) v /l + (2j ; , n &to/a )- 
sinh 2 (7r^ n i /ao) 



(45) 



An alternative way to obtain this result is o use the corresponding quantum mechanical 
reflection probability given in Now we can find the total energy released by adding 
the energy of each quantum: 



E 



1 ~ ~ . cos 2 [(7r/2)^l + (2j /)n 6t /ao) 5 

«o /=0 n = x smh {7iji :n to/a ) 



(46) 



Recall that we have obtained formulae ( ^T|) and ( f45D neglecting the terms on the right of 
equations (0), and hence these formulae give the number of particles produced due to 
the squeezing effect. In general, the acceleration effect will give additional contribution 
to the particle creation. It is important to stress again that these two types of effects give 
additive contributions to first non- vanishing order of perturbation theory. As a result the 

) as giving a part of the quanta number due to the squeezing effect 



formulae ([IT]) and ( 

are valid for more general situations, when additional contributions to the total number 
of particles have to be taken into account due to the acceleration effect. 

By using the results given above for a scalar field we can obtain the number of produced 
particles in the physically more realistic electromagnetic case. For this note that electro- 
magnetic field in a spherical cavity is a superposition of two types of modes: transverse 
electric (TE) and transverse magnetic (TM). The TE-modes correspond to the scalar 
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modes with Dirichlet boundary condition (excluding the I = mode) and TM-modes 
correspond to the scalar modes with mixed type (Robin) boundary condition 



ji(ua) + uajtiua) 



0. 



(47) 



on the sphere surface r = a. Now the number of photons produced can be presented in 
the form 



< in|iVz n |in >=< in|A^ E) |in > + < in|JV, 



r(TM) 



In 



m >, 



I = 1,2,... 



(4S 



where the first summand on the right is given by formulae fl4"T| ) and (pE5D, and the cor- 
responding expressions for the second summand are obtained from these formulae by 
replacement ji n — > ji n , where ua = ji n are solutions to the equation (E7]). 



4 Conclusion 

In this paper we have considered the particle creation from the scalar vacuum by moving 
spherical boundary. For the general case the set of equations is derived for the instanta- 
neous basis expansion coefficients and corresponding Bogoliubuv coefficients are consid- 
ered. To solve these equations and to obtain the particle number, various approximations 
can be used. Here we use the adiabatic approximation, assuming small velocities for the 
sphere motion or large frequencies for the radiation quanta. Specific examples are con- 
sidered when the number of particles produced can be explicitly found. The first one 
corresponds to the expansion or contraction of the sphere between two finite values of the 
radius, and second one corresponds to the pulse described by (f44j). The results for the 
electromagnetic field can be obtained by summing the contributions from TE and TM 
modes, corresponding to the scalar modes with Dirichlet and special Robin type bound- 
ary conditions. In general, when the acceleration effect in the particle creation cannot 
be neglected the formulae presented in previous section give only the squeezing parts of 
the total number of particles. To first non-vanishing order of perturbation theory the 
contributions from these two types of effects are additive. 
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